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Abstract 

We explore how rotational invariance is systematically recovered from calculations on hyper-cubic 
lattices through the use of smeared lattice operators that smoothly evolve into continuum opera- 
tors with definite angular momentum as the lattice-spacing is reduced. Perturbative calculations of 
the angular momentum violation associated with such operators at tree-level and at one-loop are 
presented in theory and QCD. Contributions from these operators that violate rotational in- 
variance occur at tree-level, with coefficients that are suppressed by O (a'^) in the continuum limit. 
Quantum loops do not modify this behavior in A0^, nor in QCD if the gauge- fields are smeared over 
a comparable spatial region. Consequently, the use of this type of operator should, in principle, 
allow for Lattice QCD calculations of the higher moments of the hadron structure functions. 
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I. INTRODUCTION 



Lattice quantum chromodynamics (LQCD) is a numerical technique in which Euchdean 
space correlation functions of QCD are calculated by a Monte-Carlo evaluation of the Eu- 
clidean space path integral [IJ. The computational resources are now becoming available 
for LQCD to recover the spectrum of mesons and baryons that have been observed in the 
laboratory, and to make predictions of states with exotic quantum numbers that will be 
the focus of future experimental efforts. It is also providing precise determinations of the 
matrix element of weak operators that are required to further constrain the mixing of the 
eigenstates of the weak interaction, contained in the CKM matrix. LQCD is allowing for 
a comprehensive description of the structure of nucleons, and more recently to their in- 
teractions that are crucial to the field of nuclear physics. This marks the beginnings of a 
comprehensive program to determine nuclear structure and dynamics directly from QCD. 

Space-time is pixelated, or discretized, in LQCD calculations, with the quarks residing on 
the lattice sites, and the gluon fields residing on the links between lattice sites. The lattice 
spacing, a, the distance between adjacent lattice sites, is required to be much smaller than the 
characteristic hadronic length scale of the system under study. In principle, the effects of a 
finite lattice spacing can be systematically removed by combining calculations of correlation 
functions at several lattice spacings with the low-energy effective field theory (EFT) which 
explicitly includes the discretization effects. This type of EFT is somewhat more complicated 
than its continuum counterpart as it must reproduce matrix elements of the Symanzik 
action constructed with higher dimension operators induced by the discretization [2H1]. 
While the action lacks Lorentz invariance and rotational symmetry, it is constrained by 
hyper-cubic symmetry. As computers have finite memory and performance, the lattice 
volumes are finite in all four space-time directions. Generally, periodic boundary conditions 
(BC's) are imposed on the fields in the space-directions (a three-dimensional torus), while 
(anti) periodic BC's are imposed on the (quark) gauge- fields in the time-direction. However, 
the conceptual and practical problems arising from the explicit breaking of the space-time 
symmetries of the continuum theory, down to those of a hyper-cubic lattice theory, remain 
a challenge in the continuum extrapolation of classes of observables calculated using LQCD. 
One knows, however, that as the lattice becomes finer, the full space-time symmetries of the 
continuum are in fact approximately recovered for observables involving wavelengths that 
are large compared with the scale of pixelationj^ As a result, a quantitative description of 
this restoration, as well as its implication for calculation of lattice observables, is possible. 

Efforts to reduce lattice artifacts and achieve a better behaved theory in the continuum 
limit date back to early stages of development of LQCD. Many that fall under the name 
of Symanzik improvement include a systematic modification of the action in such a way 
to eliminate O [a^) terms from physical quantities calculated with LQCD at each order 
in perturbation theory j2HH 151-115). or nonperturbatively. However, as will be discussed, 
discretization effects are known to give rise to more subtle issues; the treatment of which 
turns out to be more involved. LQCD is commonly formulated on a hyper-cubic grid, as 
a result the full (Euclidean) Lorentz symmetry group of the continuum is reduced to the 
discrete symmetry group of a hyper-cube. As the (hyper-) cubic group has only a finite 
number of irreducible representations (irreps) compared to infinite number of irreps of the 
rotational group, a given irrep of the rotational group is not irreducible under the (hyper-) 

^ For some numerical illustrations of this recovery in SU (2) lattice gauge theories, as well as the scalar 0^ 
theory, see Refs. [SHZ]- 
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cubic group. Consequently, one can not assign a well-defined angular momentum to a lattice 
state, which is generally a linear combination of infinitely many different angular momentum 
states (see for example Refs. |16nl8j ). In principle, one can identify the angular momentum 
of a corresponding continuum state in a lattice calculation from the degeneracies in the 
spectrum of states belonging to different irreps of the cubic group as the lattice spacing 
is reduced (a review of baryon spectroscopy efforts is given in Ref. [19], for some recent 
meson spectroscopy works see Refs. [201 - 123] ). However, as the density of degenerate states 
substantially increases with increasing the angular momentum, the identification of states 
with higher angular momentum becomes impossible with the current statistical precision. 
The other issue is that the cubic symmetry of the lattice allows the renormalization mixing of 
interpolating operators with lower dimensional ones. The induced coefficients of the lower- 
dimensional operators scale as inverse powers of the lattice spacing, and hence diverge as 
the lattice spacing goes to zero. Although renormalization mixing of operators is familiar 
from the continuum quantum field theory, it happens more frequently in LQCD calculations 
as the reduced symmetry of the hyper-cube is now less restrictive in preventing operators 
from mixing. To obtain useful results for, as an example, the matrix elements of operators 
from LQCD calculations, non-perturbative subtraction of the power-divergences is required 
and generally introduces large statistical uncertainties. 

To overcome these obstacles, it has been recently proposed by Dudek, et al. [2^1126] 
(and later applied to bbb system by Meinel |27| ) that by means of a novel construction of 
interpolating operators, the excited states of several mesons and baryons can be identified to 
high precision. The essence of this method is that if one uses a set of cubically invariant local 
operators which have already been subduced [28j from a rotationally invariant local operator 
with a definite angular momentum, J, while at the same time smearing the gauge and quark 
fields over the hadronic scale [291131] . the constructed operator has maximum overlap onto 
a continuum state with angular momentum J if the lattice spacing is sufficiently small. The 
subduction is assumed to be responsible for retaining "memory" of the underlying angular 
momentum of the continuum operator, while the smearing is assumed to suppress mixing 
with operators of different angular momentum by filtering contributions from ultra-violet 
(UV) modes. In another approach, states with higher angular momentum in the glueball 
spectra of 2 + 1 dimensional SU (2) gauge theories |32^ i33j are isolated by using glueball 
interpolating operators that are linear combinations of Wilson loops which are rotated by 
arbitrary angles in order to project out a particular angular momentum J in the continuum. 
In addition, the links are smeared, or blocked, in order to be smooth over physical length 
scales rather than just in the UV [34j. So by monitoring the angular content of the glueball 
wavefunction in the continuum limit with a probe with definite J, the 0~/4~ puzzle in 
the glueball spectroscopy has been tackled. The prominent feature of these works is that 
the recovery of rotational symmetry for sufficiently small lattice spacings is qualitatively 
emergent from their numerical results. 

The same issue occurs in LQCD calculations of higher moments of hadron structure 
functions, the extraction of which requires the matrix elements of local operators between 
hadronic states. Although Lorentz invariance forbids twist-2 operators with different J 
from mixing in the continuum, generally they can mix in LQCD calculations with power- 
divergent mixing coefficients [351 136] . The power-divergent mixing problem associated with 
the lower moments can be avoided by several means as described, for example, in Refs. 
|36H45j. In addition to these approaches, two methods |16| HTj have been suggested that 
highlight the idea of approaching the continuum properties of the hadronic matrix elements 
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by suppressing the contributions from the UV, and in that sense resemble the idea of operator 
smearing in the proposals described above. In LQCD calculations of non-leptonic K-decay, 
Dawson et al. ^[6] suggested that point-splitting the hadronic currents by a distance larger 
than the lattice spacing, but smaller than the QCD scale, results in an operator product 
expansion (OPE) of the currents with the coefficients of lower dimensional operators scaling 
with inverse powers of the point-splitting distance, as opposed to the inverse lattice spacing. 
This considerably reduces the numerical issues introduced by the operator mixing. In a 
different, but still physically equivalent approach, Detmold and Lin [47J showed that in the 
LQCD calculation of matrix elements of the Compton scattering tensor, the introduction of 
a fictitious, non-dynamical, heavy quark coupled to physical light quarks removes the power 
divergences of the mixing coefficients. This technique enables the extraction of matrix 
elements of higher spin twist-2 operators with a simple renormalization procedure. The 
essence of this method is that the heavy quark propagator acts as a smearing function in 
the momentum-space, suppressing contributions from the high energy modes, provided that 
its mass is much smaller than the inverse lattice spacing. 

Encouraged by the results of the numerical non-perturbative investigation of Refs. [21H26] 
and Refs. |i32| |33], as well as the results of Refs. [IHl HT], we aim to quantify the recovery of 
rotational symmetry with analytical, perturbative calculations in and QCD. In order to 
achieve this goal, we first define a composite operator on the lattice which has a well-defined 
angular momentum in the continuum limit and is smeared over a finite physical region, 
and show how the non-continuum contributions to the multipole expansion of the operator 
scales as the lattice spacing is reduced toward the continuum. Tree-level contributions to 
matrix elements that violate rotational symmetry, either by the lattice operator matching 
onto continuum operators with the "wrong" angular momentum, or matching onto contin- 
uum operators that explicitly violate rotational symmetry, scale as O (a^) as a — )■ 0. This 
includes the (naively) power divergent contributions from lower-dimension operators. In 
order to make definitive statements about the size of violations to rotational symmetry, it 
must be ensured that the tree-level scalings are not ruined by quantum fiuctuations. This is 
demonstrated by a perturbative calculation of the two-point function in scalar field the- 
ory with an insertion of such an operator. It is confirmed that quantum corrections at any 
order in perturbation theory do not alter the observed classical scalings of non-continuum 
contributions. This result is comparable with finite size scaling results of the leading irrele- 
vant operator that breaks rotational invariance in three dimensional O (N) models given in 
Refs. |18| H9] . The critical exponent p introduced there, has a realization in terms of small-a 
scaling of the leading rotational invariance violating terms in this calculation. Its value is 
shown to be consistent with the results presented here. 

After gaining experience with this operator in scalar field theory, the generalization to 
gauge theories is straightforward. Special attention must be paid to the gauge links that 
appear in the definition of gauge-invariant operator(s) that are the analogue of those con- 
sidered in the scalar field theory. Also, it is well known that the perturbative expansion of 
operators used in LQCD are not well-behaved due to the presence of tadpole diagrams [50) . 
Naively, tadpoles make enhanced contributions to the matrix elements of the operators we 
consider, and that tadpole improvement of the gauge links and smearing of the gluon fields 
are crucial to the suppression of violations of rotational symmetry. After discussing the 
continuum behavior of the QCD operator (s), and their potential mixings, which violate ro- 
tational invariance at O(a^), we determine the renormalization of the operator (s) on the 
lattice at one-loop order. The leading rotational invariance violating contributions to the 
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FIG. 1. A contribution to the lattice operator defined in Eq. ([T]), with |n| < A^. All the points 
inside the three-dimensional spherical shell |na| = Na are included in the operator. The two length 
scales defining the operator, the lattice spacing, a, and the operator size, Na = 1/A, are shown. 

renormalized lattice operator are suppressed by 0{asa'^), (where ag = g1/{'iTi) and Qs is the 
strong coupling constant) provided that the gauge fields are also smeared over a physical 
region similar to the matter fields. This means that the leading rotational invariance vio- 
lating operators introduced by the quantum loops make subleading contributions compared 
to tree-level, 0{a?). The loop contributions that scale as 0{asa) do not violate rotational 
symmetry, and hence are absorbed into the operator Z-factor. 



II. OPERATORS IN SCALAR FIELD THEORY 

The goal is to construct a bilinear operator of the scalar fields on a cubic lattice which has 
certain properties. First of all, as it was discussed earlier, it has to be smeared over a finite 
region of space. This physical region should be large compared to the lattice spacing, and, 
for our purposes, small compared to typical length scale of the system to allow for a pertur- 
bative analysis. The spatial extent of the operator can be identified with its renormalization 
scale. Secondly, it is required to transform as a spherical tensor with well-defined angular 
momentum in the continuum limit. An operator that satisfies these conditions is[^ 

3 '"'-^ 

eL,M (x; a, N) = ^ </. (x) </) (x + na) Yl,m (n) , (l) 

n 

where n denotes a triplet of integers, and it is normalized by the spatial volume of the region 
over which it is distributed. 0(x) is the scaler field operator, is the maximum number 
of lattice sites in the radial direction, and Y^^j (fi) is a spherical harmonic evaluated at the 

^ This corresponds to one particular choice of radial structure of the operator. However, the results of the 
calculations and the physics conclusions presented in this work do not change qualititively when other 
smooth radial structures are employed, such as a Gaussian or exponential. 
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angles defined by the unit vector in the direction of n, n, as shown in fig. [T| This operator 
can also be written in a multipole expansion about its center as 



9l,m (x; a, N) = — ^ (x) (an ■ V)'= (x) Yl,m (A) , (2) 

n k 

where the gradient operator acts on the x variable, V = Vx- 

Although the operator O^^m (x; a, N) is labeled by its angular momentum in the contin- 
uum limit, from the right hand side of eq. (|2|, it is clear that it is a linear combination 
of an infinite number of operators with angular momentum compatible with its parity. To 
be more specific, consider the M = component of the operator expanded in a derivative 
operator basis, 

e..(x;a.iV) = ^£%<^ol?,(x;a) , (3) 

L',d 

where O^^, (x; a) are defined in Appendix A The operator subscript denotes that there are 
L' free indices in the derivative operator, wmle d denotes the total number of derivatives. As 
is discussed in the Appendix |A| there are operators in this basis which are not rotationally 
invariant but only cubically invariant. C^q.j^,q (N) are coefficients of each operator in the 
expansion whose values are determined by matching eq.([2]) with eq. ([s]). Finally A = l/(Na) 
is the momentum-scale of the smeared operator which is kept fixed as the lattice spacing 
is varied. Therefore, as the lattice spacing decreases, more point shells (shells of integer 
triplets) are included in the sum in eq. The convergence of this derivative expansion is 
guaranteed as the scale A is set to be much larger than the typical momentum encountered 
by the operator. 



A. Classical Scalar Field Theory 

In order for the operator to recover its continuum limit as the lattice spacing vanishes, 
the coefficients C^^q.j^,q should have certain properties. First of all, those associated with the 
operators with L ^ V as well as the rotational invariance violating operators, should vanish 
as a — 7- 0. Also the coefficients of rotational invariant operators with L = L' should reach 
a finite value in this limit. These properties will be shown to be the case in a formal way 
shortly, but in order to get a general idea of the classical scaling of the operators and the 
size of mixing coefficients, we first work out a particular example. Consider the operator 
Osfi (x; a, N) expanded out up to five derivative operators, 

^3,„(x;a,A') = %i^O<"(x;a) + %5^5o?'(x;a) + !%^5o?'(x;a) + 

(x; a) + H%^0<£ (x; a) + %™og) (x; a) + (4) 

%™OSL(x;a) + 0(S) , 

where the superscript RV denotes the rotational invariance violating operator and its corre- 
sponding coefficient in the above expansion. 
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FIG. 2. The tree-level values of the coefficients c'^.^i^ appearmg in eq. (4) as a function of the 
largest n-shell included in the summation in eq. ([T]). 



The numerical values of the coefficients in eq. Q, at the classical level, as a function 
of the maximum shell included in the sum in eq. ^ are shown in fig. [2] and fig. [s] From 

these plots it is clear that while the coefficients Cfj^^ao and Cf^.^Q reach a finite value for 
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N 

FIG. 3. A comparison between the tree-level coefficients C^'q-i^/q to illustrate the relative rates of 
convergence to the continuum limit. 

large N, the coefficients of lower and higher angular momentum operators, as well as the 
rotational invariance violating operator, approach zero. To find the values of the leading 
order (LO) coefficients in this limit, as well as to see how the non-leading contributions scale 
with = l/(Aa), one can apply the Poisson re-summation formula to the right hand side 
of eq. g. 



9,,M (x; a, AT) = ^ ^ ^ E / ^ " ^) <^ (^) " <^ (^) ' 

fc ■ p 

(5) 

where p is another triplet of integers, and the p summation is unbounded. The continuum 
values of the coefficients obtained in the N ^ oo limit, corresponding to the p = term in 
eq. ([5]), are 



15 7 £-1 



4 V vr + 



with (i = 3, 5, ... , 



(6) 



while the other coefficients in eq. Q vanish in this limit as expected. The LO corrections 

to these continuum values can be calculated as following. The deviation of C30.30 from its 
continuum value can be found from 



'30 



3 {No) 
4^ 3! 



J2 f'dyy'dnye^'-'^^-y<P{^) (y-Vf0(x) ^3,0 (y) 



(7) 



where V = V^Cz and the variable in eq. ([T]) is redefined to lie between and 1, and it is 
straightforward to show that 



6C. 



(3) 
30:30 



1 



A^2 327r2 V TT 



p^O 



cos(27rA^|p|) / 3 



I u I T r I 1 2 4 6 

Pl + 15 IpI - —p^ 



25 

y 



(8) 
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It is interesting to note that, after trading for l/(aA), the finite lattice spacing corrections 
are not monotonic in a, but exhibit oscillatory behavior, which is clearly evident in fig. [2j 

The deviation of Cgp^o from its continuum value of zero follows similarly, and is found to 
scale as ~ 1/^^, 



5C. 



1 3 /7 ^ cos(27rAr|p|) , 4 4. 



= 

30;10 - ^2 15^2 y ^ 



As in the case of the operator that conserves angular momentum in the continuum limit, the 
sub-leading correction (and in this case the first non-zero contribution) to the coefficient is 
suppressed by l/N"^. This can be shown to be the case for all the sub-leading contributions 

to the coefficients C^jvf L'M' follows. As is evident from eq. (5), the integrals that are 
required in calculating deviations from the continuum values have the general form 



(10) 



which can be written as 



Tl\...lk 



3 {Naf 1 



4vr k\ {i2'KN) 
3 (Na)^ Am^ 



k 



d 


d 


dPh 


'dPi, 


d 


d 


dpi,' 


'dPi, 



cos(27rAf|p|) sm(27rAf|p|) 
{2nN\p\f {27rJV|p|)2 

higher orders in 1/A^, depending on whether L is even or odd. Thus the LO contribution from 



The y integration over the Bessel function gives rise to either — ™! Z, fJ or — '^'^ Z, ; up to 

° ° (27rAf|p|)^ {27rAf|p|)^ ' ^ 



eq. (11) in the large limit is obtained by acting on the numerator with the p derivatives, 
producing k powers of A^, multiplying the 1/A^^ from the denominator. Therefore, eq. (11) 
scales as 

1 1 1 

pi-ik ^ (Nat—- -— , (12) 

and, in general, the deviation of any coefficient from its continuum value is suppressed by 
1/A^^ = A^a^. This result implies that in calculating the matrix element of L = 3 operator, 
one has a derivative expansion of the form 

A3^3,o (x; a, AT) = ^0«(x;a) + a, (x;a) + «3 (x; a) 

+ «4 ^OP^)(x;a) + a5 0E(x;a) + a, (x; a) 

+ a4^^-- ^) + ^ (5) ' (1^) 

where the mixing with L 7^ 3 operators (with coefficients 01,2,3,7,...), as well as the operator 
with broken rotational symmetry (with coefficient a^), vanish in the large A^ limit, while 
the coefficients of L = 3 operators (with coefficients a5,6,...), are fixed by the scale of the 
operator, A. It is clear that for A^ = 1 and A = 1/a, where no smearing is performed, the 
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problem with divergent coefficients of the lower dimensional operators is obvious, as, for 
example, the coefficient of O^P (x; a) diverges as as a — )■ 0, as is well known. 

The fact that all the sub-leading contributions to the classical operator are suppressed 
at least by regardless of L and V can be understood as follows. In the classical limit, 

where the short distance fluctuations of the operator are negligible, the operator does not 
probe the distances of the order of lattice spacing when a — )■ 0. The angular resolution of 
the operator is dictated by the solid angle discretization of the physical region over which 
the operator is smeared, and therefore is proportional to l/N"^. The question to answer is 
whether the quantum fluctuations modify this general result. 

Before proceeding with the quantum loop calculations, it is advantageous to transform 
the operator into momentum-space to simplify loop integrals. This can be done easily by 
noting that for zero momentum insertion, the operator acting on the field with momentum 
k is 

hM{k;a,N) = —— e^''^^ Y^m (n) (k) 0(-k) , (14) 



which, after using the partial-wave expansion of e*'' "'^ and the exponential term resulting 



from the Poisson relation, can be written as 



(k; a, iV) = 6v^ (k) (-k) E \ ^'^^ + f ^ + 



2L + 1 

p Li,Mi,L2,M2 

X (LiO; L2O \L0) (LiMi; L2M2 \LM)Yl,m, i^j^) Yl,m, (f^p) 

X f rfy JL, (aAr|k|y) ji, (27riV|p||/) . (15) 

Although this form seems to be somewhat more complicated than in position-space, it turns 
out that it is advantageous to work in momentum-space when dealing with higher angular 
momenta, as well as for M 7^ 0. Further, the dimensionless parameters |k|/A and A^ that 
define the physics of such systems are now explicit. It is straightforward to show this form 
recovers the values of the leading and sub-leading coefficients given in eq. (|8| and eq. (|9|), 



and it is worth mentioning how they emerge from eq. (15). For a non-zero value of |p| and 
N = 00, the spherical Bessel function j^^ {2ttN |p| y) vanishes for any value of L2. However, 
for large values of A^ but |p| =0 the only non-zero contribution is from L2 = 0, and 
thus Li = L, leaving a straightforward integration over a single spherical Bessel function 
ji (aN \\i\y) to obtain the continuum limit given in eq. (|6|. Extracting the subleading 
contributions and the violations of rotational symmetry is somewhat more involved, and we 
provide an explicit example in Appendix [B} 



B. Quantum Corrections in Xc/)^ 

In order to determine the impact of quantum fluctuations on the matrix elements of 6l^m, 
defined in eq. ([T|, we consider loop contributions in theory. Beside its simplicity which 
enables us to develop tools in performing the analogous calculations in Lattice QCD, this 
theory corresponds to some interesting condensed matter systems. For example, three di- 
mensional 0(N) models, which describe important critical phenomena in nature, have a 
corresponding field theory formulation. As pointed out in Refs. [HI HQ], anisotropy in 



11 



FIG. 4. One-loop correction to the two-point function with an insertion of 6l,m in ^fp^ 



space either due to the symmetries of the physical system, or due to an underlying lattice 
formulation, will result in the presence of irrelevant operators in the effective Hamiltonian 
which are not rotationally invariant, and introduce deviations of two-point functions from 
their rotationally invariant scaling law near the fixed point. However, as the rotationally 
invariant fixed point of the theory is approached, the anisotropic deviations vanish like 1/^^ 
where is the second moment correlation length derived from the two-point function, and 
p is a critical exponent which is related to the critical effective dimension of the leading 
irrelevant operator breaking rotational invariance. It has been shown that in the large N 
approximation of O {N) models, p ~ 2 for cubic-like lattices. In the following, it will be 
shown that, by inserting 9l,m defined in eq. ([l]) into the two-point function, the same scaling 
law emerges when approaching the rotational-invariant continuum limit of theory. 
At tree level, the contributions to the two-point function from an insertion of 6l,m at 



zero momentum transfer has been already discussed in section II A At one-loop, there is 
only one diagram with an insertion of 6l,m that contributes to the two-point function, as 
shown in fig. |4j This diagram introduces corrections only to the L = matrix element as 
there are no free indices associated with the loop. The lattice integral associated with this 
one-loop diagram is 



J, 



LM 



3A 

AttN^ 



|n|<Af 



jkna 



(16) 



where k"^ 



2 



A is the coupling constant and m is the mass. The three- 



momentum integration can be evaluated by noting that the region of integration can be 
split into two parts: region I where < |k| < vr/a and therefore is rotationally symmetric, 
and region II where n/a < |k| < y/Sir/a which consists of disconnected angular parts. Also 
as the three-momentum integration is UV convergent, a small a expansion of the integrand 
can be performed. Using eq. (15), the contribution from region I to the p = term in the 
Poisson sum is 



4Z (p = 0) 



3A 
(2^ 



3A 
16^ 



dkA 



dkP / da 



1 



k X 2 

f^2 _j_ ^2 



dy y'^ Jl {aN |k| y) 



(17) 
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where 



Jlm — 2V7r(^L,o(^M,o 



/Aa r 1 



dqq^ 



[q^ + ql + m^/K'Y Jo 



f 

Jo 



dy jo (qy) 



J 



NLO 
LM 



\ I Aa I Aa 



X 



Sl.o S 



L,0 (JM,0 



[q^ + qj + myA^Y 
■1 

dyy^jo (qy) 





dy y"^ j4 {qy) 



,(18) 



with q = |k|/A and q^ = k^/A. The LO integral, Jf^, is convergent, while the NLO 
contribution, J^^'^, while not convergent, is not divergent, but is of the form sin (Nir) /N'^. 
This implies that they depend on the ratio of the two mass scales, A and m, but without 
inverse powers of a. So as a — > 0, the LO L = operator makes an unsuppressed contribution 
to the L — matrix element, while the contributions to this matrix element from the NLO 
rotational-symmetry violating L = and L = 4 operators are suppressed by 1/N'^. 

A simple argument shows that contributions from integration region II, for which n/a < 
|k| < y/Sir/a, are also suppressed by l/N"^. After defining a new momentum variable 
I IX — k^a and l'^ — + + 1^, the p = term of the Poisson sum in region II i 



IS 



Ji'S (P = 0) 



3A 



7r r 
dl4 I 

■TT J TT 



\/37r n 

dl /2 / dVti 

Ylm (^i) 
4E;.sin2 {lj2) + a?m^)' 



f'dyy'hiNly) , (19) 

^0 



where / (fli) identifies the angular region of integration, and whose parametric form does 
not matter for this discussion. This region still exhibits cubic symmetry, and gives rise to 
contribution to the L — 0,4,6,8,... operators. On the other hand, the three-momentum 
integration is entirely located in the UV as a — > 0, and thus 



sin^ (/i/2) + sin^ {I2/2) + sin^ {1^/2) + sin^ {U/2) > 1 



(20) 



Also, integration over the Bessel function brings in a factor of — cos (Nl) /(N'^l^), up to 
higher orders in So the integrand does not have any singularities in region II of the 

integration, and is bounded. As a result. 



(21) 



and consequently J^l^ (p = 0) itself is suppressed by This completes the discussion 

of the p = term in the Poisson sum, corresponding to a zero- momentum insertion of the 
continuum operator into the loop diagram. It then remains to determine the scaling of the 
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p 7^ terms in the summation in the large hmit. The integral arising from the p 7^ 
terms is, up to numerical factors, 



„ /■« d^fc 1 
Xp^o ~ / . / T^TTi-. r2 ^L^Mi i^k) YL2M2 i^p) 

X [ dyy' jL, {Na |k| y) jl, {2'kN |p| y) . (22) 
^0 

This integral is finite in UV, and integrand can be expanded in powers of a, giving a leading 
contribution of 

t^o-i-f-a (2^) (g^ + gi + "^VA2) 

X f dyy^ jL^qy) JLA'^'^N\p\y) . (23) 
Jo 

A non-zero angular integration requires that Li = 0, and the integral is suppressed at 
least by a factor of l/N"^ as integration over the Bessel functions introduces a factor of 
1/ {27tN IpD^ up to a numerical coefficient and a bounded trigonometric function at leading 
order in The next order term in the small a expansion of the integrand can be easily 

shown to bring in an additional factor of l/N"^. So one can see that the p 7^ terms in 
the Poisson summation, which give rise to non-continuum contributions to the two-point 
function at one loop, are always suppressed by at least a factor of l/N"^. 

The result of the one-loop calculation is promising: all the sub-leading contributions that 
break rotational symmetry are suppressed by l/N"^ compared to the leading L = continuum 
operator contribution to the two-point function. A little investigation shows that this scaling 
also holds to higher orders in A0'^ theory. Suppose that the operator is inserted into a 
propagator inside an n-loop diagram contributing to the two-point function. Considering 
the continuum part of the operator first, the leading term in the small a expansion of the 
integrand gives rise to 2n propagators, while the integration measure contributes 4n powers 
of momentum. Although this appears to be logarithmically divergent, the spherical Bessel 
function contributes a factor of inverse three-momentum and either a sine or cosine of the 
three-momentum, rendering the diagram finite. The same argument applies to the NLO term 
in the small a expansion of the integrand, resulting in a suppression of the breaking 

of rotational invariance. Insertion of the non-continuum operator in loop diagrams are also 
suppressed by l/N"^ for similar reasons. 

The interpretation of finite size scaling results presented in Refs. |18| HQ] in terms of 
what has been observed in this section is now straightforward. Near the critical point, the 
correlation length is the only relevant physical scale in the problem, and tends to infinity. 
So as the critical point is approached, one does not probe the underlying lattice structure 
as the correlation length becomes much larger than the lattice spacing, and extends over 
an increasing number of point shells. In comparison, inserting an operator which only 
probes distances of the order of a physical scale that is much larger than the lattice spacing, 
resembles the physics near a rotational-invariant fixed point, and the same scaling law for 
the non-rotational invariant operators is expected (in the same theory) as the lattice spacing 
goes to zero. 
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FIG. 5. One-loop QCD corrections to the fermionic two-point function with an insertion of 6l. 



M, 



given in eq. (24), at zero external momentum 



III. OPERATORS IN QCD 



The necessity of introducing a gauge link to connect the fermionic fields in a gauge invariant 
way, makes the discussion of the operator and its renormalization more involved in gauge 
theories. The reason is two-folded: firstly as is well known, perturbative LQCD is ill- 
behaved as a result of non-vanishing tadpoles which diverge in the UV, making the small 
coupling series expansion of the operators slowly convergent. The other difficulty is that 
as the operator is smeared over many lattice sites, the links are necessarily extended links. 
Thus, to analytically investigate the deviations from a rotational invariant path, working 
with a well-defined path on the grid is crucial. In this section, the strategies to deal with 
these problems are discussed, and the scaling laws of different operator contributions to the 
two-point function in QCD with an insertion of the smeared operator are deduced. 

In position-space, perhaps the simplest gauge-invariant smeared operator of quark bilin- 



ears is 



|n|<7V 



9l,m (x; a, N) 



ip (x) t/ (x, X + na) ip {■X + na) Yl^m (n) 



(24) 



where 



[/ (x, X + na) = e 



/■x+na 

igj:+'^^A{zydz ^ l+ig A (z) ■ dz + O (q^) 

J X 



(25) 



where the actual path defining U will be considered subsequently. As the fermion operator 
is a spin singlet, S = 0, the total angular momentum of this operator in the continuum is 
J = L. One could also consider operators of the form 



0JL,M (x;a,iV) 



4vrA^3 



|n|<Af 



x) ?7 (x, X + na) -i/^ (x + na) Yl^m (n) 



(26) 



which can be used to form operators with J = L + 1, L, L — 1. It is clear that the set of 
operators with angular momentum J will mix under renormalization, but the vector nature 
of QCD precludes mixing between the ipip and operators in the chiral limit. However 

to capture the main features of operator mixing in the continuum limit of LQCD, it suffices 
to work with the simplest operator, in eq. (24). At tree-level, the contributions of this 



operator away from the continuum limit scale in the same way as in the scalar theory, with 
contributions that violate rotational invariance suppressed by ~ l/N"^. 

Let us first discuss the one-loop renormalization of the operator in the continuum. There 
are three one-loop diagrams contributing to the operator renormalization as shown in fig. |5| 
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The diagram in fig. |5^ results from inserting tlie leading order term in the small coupling 
expansion of the operator in the loop. At zero external momentum this diagram is 

r(-) ~ -T'^T^^^ f dyy^ [ dQ, [ J^^ ^^i^f^^r + f l\ .....yy ^^T) 

which is clearly convergent in the UV. Also it contains L = as well as L = 1 operator as 
can be seen from the angular part of the integral 

^ f dQydQ^[f,{e,m,k,)+f2{e,m,k,)k-^] YL>M'{^^.) >?'m' (^y) Y^m i^y) 

L',M' 

47r/i [k^, k4,m) 5l,o5m,o 



in 

+ Y {k^,k4,m) \k\SL, 



7i 7= +^72 7= + Ishifl 



V2 ) "\ V2 



(28) 



where fi and /2 are some functions of their arguments. One can check however that as 
m/A — )■ (the chiral limit), the contribution to the L = 1 operator is suppressed by the 
quark mass. 

The diagrams in fig. ^jp comes from the next term in the expansion of eq. (25). It is 
straightforward to show that the Feynman rule for the one-gluon vertex with zero momentum 
insertion into the operator is 

^/ = ^ E aan' (iT-poTW (^''''''"" - ^"'■"") ^^'^ ' 

where the radial path between points x and x + na is taken in evaluating the link integral, 
p and p' are the momenta of incoming and outgoing fermions respectively, A is the Lorentz- 
index of the gluon field, and k is the momentum of the gluon coming out of the vertex. 
Note that in principle, any path between points x and x + na can be taken in the above 
calculation, but if one is interested in deviations of the renormalized lattice operator from 
the rotational invariance compared to the continuum operator, a path between two points 
should be chosen in the continuum in such a way that it respects rotational invariance 
explicitly. Any path other than the radial path, on the other hand, is equivalent to infinitely 
many other paths resulting from rotated versions of the original path around the radial 
path. To reveal rotational invariance at the level of the continuum operator, an averaging 
over these infinite copies of the path is needed, and this makes the calculation of the link 
more involved. 

Now at zero external momentum, using expression (29) with p = 0, the contribution from 
the second and third diagrams in fig. [Sjo is 



(27r)^ (F + m2)A;2 



X ^ (e^^^'^-y - 1) Yr^M {^y) ■ (30) 

As is evident, because of a non-oscillatory contribution to the operator, there is a loga- 
rithmically divergent piece from the above integration contributing to the L = operator, 
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a 



b 



FIG. 6. The tadpole contribution consists of the conventional tadpole diagram (a), which vanishes 
when using a mass-independent regulator in the continuum (such as dimensional regularization) , 
as well as the diagram shown in (b) which is of the order of a^/ | Axj'^, where Ax is the distance 
between two gluon vertices. 



which along with the logarithmic divergent contribution from wavefunction renormahzation, 
contributes to the anomalous dimension of the operator. Also the angular integration of the 
above expression: 



/ 



dflydflk 



1 , y -7 

1 + — m 

iky 



J dily [g, {Nay |k|) + (Nay |k|) my • 7] Y^m i^y) , (31) 



indicates that as before, in addition to L = operator, an L = 1 contribution is present 
which is finite at UV, and can be shown to vanish for m/A — )■ 0. gi and g2 are some functions 
of their arguments whose explicit forms do not matter for this discussion. 

The last diagram in fig. [s] corresponds to the O {g'^) term in the small coupling expansion 
of the gauge link. It contains the tadpole of the continuum theory whose value depends 
in general on the regularization scheme. For example, by using a hard momentum cutoff 
which is matched easily with the lattice regularization, it diverges quadratically. However, 
it is not hard to see that in dimensional regularization which respects the full rotational 
symmetry of the continuum, it vanishes in d = 4, therefore it does not contribute to the 
renormahzation of the continuum operator. But the fourth diagram in fig. [5] does not only 
include the conventional tadpoles, fig. |6^, it also contains the diagram where a gluon is 
emitted by the Wilson line inside the operator and then absorbed at another point on the 
Wilson line, fig. |6]3 as a consequence of the matter fields being separated by a distance na. 
It is straightforward to show this diagram is convergent, and scales by a^/ |Ax|^ where Ax 
is the distance between two gluon vertices and is evaluated at the energy scale of the 
order of 1/ |Ax|. This completes the qualitative discussion of the operator renormahzation 
and mixing at one-loop order in the continuum. 

Let us start the discussion of the lattice operator by assuming that its definition is still 



given by eq. (24). However, this can be shown to be a naive definition of the operator on the 
lattice. The reason is implicit in the discussion of tadpoles given above. Although tadpoles 
are absent from the operator renormahzation in the continuum, on the lattice, they are 
non- vanishing, and result in large renormalizations, as can be seen in perturbative lattice 
QCD calculations. As was suggested long ago by Lepage and Mackenzie [50], to make the 
perturbative expansion of the lattice quantities well-behaved, and to define an appropriate 
connection between the lattice operators and their continuum counterparts, one can remove 
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Na 



FIG. 7. Tadpole diagrams contributing to the smeared operator at one-loop order. Shown in the 
right are the number of diagrams of each type. 

tadpoles from the expansion of the lattice operators in a non-perturbative manner by dividing 
the gauge link by its expectation value in a smooth gauge, 

U {x,x + afi) — 7- — U {x, X + ajl) , (32) 

Mo 

where a simpler, gauge invariant choice of Mq uses the measured value of the plaquette 
in the simulation, Mo = (|Tr (t/p^a^))^^^. There remains still another issue regarding the 
tadpole contributions to the smeared operator which is not fully taken care of by the simple 



single-link improvement procedure explained above. The operator introduced in eq. (24) is 
smeared over several lattice sites, and as a result includes extended links. As will be explained 
shortly, in spite of O {as) corrections due to tadpoles from a single link, there is an O [Nas) 
enhancement due to the tadpoles from the extended link with length ~ Na. So although a 
non-perturbative tadpole improvement could introduce non-negligible statistical errors, this 
improvement is crucial, otherwise the relation between the lattice smeared operator and the 
corresponding continuum operator is somewhat obscure. 

The reason for the O (Nag) enhancement of tadpoles from the extended links can be 
illustrated by working out a particular example. Suppose that the link is extended between 
points X and x -|- Naii entirely along the 1 axis. Then in order to make a tadpole, not 
only can each gauge field be contracted with the other gauge field belonging to the same 
elementary link, but also it can be contracted with a gauge field from one of the remaining 
A^ — 1 elementary links (see fig. [T]). Note that each diagram in fig. [T] comes with a multiplicity 
of N — m, where m is the number of links between the contracted gluonic vertices. At LO 
in a, the corresponding contribution from the extended tadpole (ET) is of the form 



r(-) ~ a.a- I " d^k , ' , , ~ S , (33) 



fc? ~l~ Ajn -|- k'i -\- k'i 



m 



2 



from which the contribution from all the diagrams in fig. [7] can be obtained. 



N-1 

V(Ar-m)^ = 0{Nas) . (34) 

m=l 
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Note that the m = term, corresponding to the first diagram in fig [7} has been excluded 
from the above sum as it is just the single link tadpole contribution. Given that there are 
N single links, the total contribution from single link tadpoles is O (Nag) as well. 

Another issue with the extended links is the fact that without tadpole improvement, 
breakdown of rotational symmetry occurs at O {Nag). The reason is that without tadpole 
improvement of the extended links, contributions from the different Ai irreps in a given 
point shell are normalized differently. For example, there are more tadpole diagrams at 
O {g"^) contributing to an extended link between points (0, 0, 0) and (2, 2, 1) (six single links) 
than to an extended link between points (0, 0, 0) and (3, 0, 0) (three single links) although 
both points belong to the same point shell (i.e. have the same separation in position space). 
This fact magnifies the necessity of tadpole improvement as well as providing a prescription 
for an appropriate improvement of an extended link. As the expectation value of a link 
belonging to a given Ai irrep in a given shell is in general different from the expectation 
value of the link belonging to another Ai irrep in the same shell, one needs to redefine the 
link in a given irrep by dividing it by its expectation value in the same irrep. 



llAt (x, X + an) — 7- 



-Uy^i (x, X + an) 



(35) 



where ti^i = \Ua{ {x,x + an)j, and the A^s are different Ai irreps belonging to the n^- 
shell. With this prescription for tadpole improvement of the extended links, the renormalized 
operator is assured to be safe from large rotational invariance breaking effects of the order 
of O {Nag). With this new definition of the gauge link, eq. (24) is now a well-defined lattice 
operator with an appropriate continuum limit which can be used in our subsequent analysis. 

As the cancellation of the tadpole diagram is assured by the new definition of the operator, 
there are only three one-loop diagrams that contribute to the renormalization of the lattice 
operator. The first diagram in fig. [5] corresponds to the following loop integral at zero 
external momentum for Wilson fermions, 



(2vr) 



7pCos 



kpQj 



ir sm 



kpQj 



-^E,7.^ + M(fc) 



X 



X 



■y^ cos — ir sm 

' ' 2 / V 2 



;?E.sin^(¥) 



^ YlM i^r 



(36) 



where M (k) = M + 2r/a^sin^ {k^a/2), and r is the Wilson parameter. Clearly at LO in 

the lattice spacing, one recovers the corresponding diagram with the insertion of the contin- 
uum operator, eq. (27), and so it contributes to both the L = and L = 1 operators. Note 
that although the integration region is not rotationally symmetric like the continuum inte- 
gral, the convergence of integral at UV ensures that the contributions from non-rotationally 
symmetric integration region II, defined in section |IIB[ are suppressed by additional powers 
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of compared to the rotational invariant region I: 



167r4 



did 



dl V 



dyy^jL (Nly) 



(37) 



where: = k^a and = Z^ + Z^ + Zg. The integrand is clearly convergent, and the integration 
region is entirely in the UV, and so the only dependence on a = 1/(AA^) comes from 
the integration over the Bessel function, giving a LO contribution proportional to l/N"^. 
However, the first sub-leading contribution from this diagram scales as ~ ag/N for Wilson 
fermions instead of ~ as/N"^- The reason is that the small a expansion of the integrand 
in eq. (36) includes terms at O (a) which is proportional to the Wilson parameter. The 
integrand scales as ~ 1/k^ multiplied by the spherical Bessel function in the UV which still 
gives rise to a convergent four-momentum integration for any value of L, 



^■p(5a,r) 



~ a 



d^k 



k^ 



dy 3l {Naky) 



aA 



1 

N 



(38) 



These contributions are rotational invariant, and will be included in the renormalization Z- 
factor of the operator when matching the lattice operator with its continuum counterpart. 



Further, the integrals that appear at O (a^) in an expansion of eq. (36) are also convergent, 



and the terms containing rotational invariance breaking contributions are suppressed by 
l/N"^ . This completes discussion of the first one-loop diagram of fig. [sj 

The second diagram contains the one-gluon vertex operator, and requires evaluating a 
line integral over the path on the grid defining the extended link. As was pointed out in the 
discussion of the path in the continuum, in general any path can be chosen in evaluating the 
operator both in the continuum or on the lattice, but requiring the recovery of rotational 
symmetry at the level of the operator means that the extended link has to exhibit rotational 
symmetry in the continuum limit. As already discussed, the simplest rotational invariant 
path in the continuum is the radial path between the points, so it makes sense to try to 
construct a path on the grid which remains as close as possible to the radial path between 
points X and x + na as it passes through the lattice sites. One might expect though that 
choosing a path in continuum which is the same as its lattice counterpart is a more legitimate 
choice. One example of such a path is an L-shaped path. However, it is not hard to verify 
that the L-shaped link does not restore rotational invariance in the continuum limit as the 
continuum path explicitly breaks rotational symmetry. So the problem of evaluating the one- 
gluon vertex of the smeared operator is reduced to finding the closest path to the straight line 
on the grid. In a lattice calculation, one can, in principle, construct an algorithm which finds 
a path on the three-dimensional grid in such a way that the area between the path and the 
rotational invariant radial path is a minimum. One such algorithm has already been used in 
Ref. [32j to construct a path that follows the straight line between sites A and B as closely 
as possible, by forming a diagonal link at each step which has the maximum projection 
onto the vector A^. By this construction of "super-links", the authors have been able to 
form arbitrary (approximate) rotations of the Wilson loops, therefore constructing glueball 
operators which project onto a definite spin J in the continuum limit. However, the analytic 
form of the super-link has not been given. In appendix [Dl a method to evaluate the link 



on such a path is illustrated with a small number of examples. For the following discussion 
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x+na x+na x+na 




X x X 

a b c 

FIG. 8. a) The link between points x and x + na for n = (2, 1, 0) which remains as close as possible 
to the diagonal link, b) The link between the same points for n = 2 (2, 1,0) which consists of two 
separate links of part a) with the lattice spacing being halved, c) The link for n = 2^ (2, 1, 0) which 
consists of 2^ separate links of part a) with the lattice spacing divided by 2^ . 



however, a particular example has been considered which encapsulates the essential features 
of the recovery of the rotational path, and gives us an idea how to deal with the general 
case. 

Suppose that the link connects points x and a; + na on a cubic lattice where n = 
^ (Q,1,0), and ao = 2^a. As usual a denotes the lattice spacing, and Q is an arbitrary 
integer. The continuum limit is recovered when the integer K tends to infinity for a finite 
value of ao- Then as is shown in appendix |D| for a path which is symmetric under reflection 
about its midpoint and remains as close as possible to the vector na (see fig. [S]), the O (g) 
term in the momentum-space expansion of the link has the following form 



(g) 



^ nK 



Sin {^) 



sm 



( q na \ 



(q) 



[Ay (q) 

sin (gg,.ao/2 ^+^) 
sin(ga;ao/2 



K+l] 



COS 



2^+2 ~^ 2^~^^ 



(39) 



As — 7- oo limit which corresponds to a — t- 0, one obtains 



f/(^f) (g) = 2i^e^'i-"''/2 



sm 



(^) 



q ■ na 



A ■ na + — {q^Q + qy) A na 



a 

"24 



QA^ao [ql _ 1) + ^Qq,^qy + 3^2) + q (^4^) 



(40) 



recovering the continuum link, given in eq. ( |29[ ), and contains broken rotational invariance 
contributions which are suppressed by ~ 0{a^. This scaling has been shown in appendix 
|D|to hold for vectors n of the forms: f (Q, 1, 1), f (Q, Q, 1) and f (Q, Q, Q) as well. 

Let us now examine how the insertion of this contribution from the operator modifies the 
scaling of the rotational invariance violating operators at one-loop. The contribution from 
the second diagram in fig. [5] with the insertion of this vertex can be calculated order by 
order in small a by expanding the vertices and propagators as before. At the LO one gets 



p(5fe) 



- 1 



47riV3 (27r)^ {k^ + m^) k^ tk ■ na ^^"^ 



an ■ 7 + — {k^Q + kyf an • 7 - —Q [kl (Q^ - l) + SQk^ky + 3A;J) -f^ao 



24 



.(41) 
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Clearly, after adding the contribution from the third diagram in fig. |5| the LO contribution 



from the above expression, the first term in the bracket of eq. (41), recovers the results ob- 
tained previously for the insertion of the continuum operator, up to suppressed contributions 
from the integration region II, as discussed before. Therefore this term contributes to the 
L = operator with a logarithmically divergent coefficient, which along with the wavefunc- 
tion renormalization contributes to the anomalous dimension of the lattice operator. Note 
that the wavefunction renormalization gives rise to a logarithmically divergent contribution 
to the L = operator at LO in the lattice spacing, recovering the continuum result, and 
the sub-leading contributions are suppressed at least by a = 1/(A^A) for Wilson fermions. 
This term also contains and L = 1 operator which is proportional to m, and vanishes in the 
chiral limit. 

The second term in the bracket of eq. ( [4T| is O(a^), and can be written as 

-|- 77^2) 

X {k,Q + kyf Ylm {^n) 

- O (/a°) . (42) 

This scaling arises as a result of the UV divergence of the non-oscillatory contribution to 
the integral and is entirely a UV effect. For this term there is no dependence upon n and as 
such the factor of is canceled by a corresponding A^^ from the sum. Terms proportional 
to the mass are convergent in the UV, and as such are suppressed by in the continuum 
limit. 



The last term in the above expression eq. (41 ) contains rotational breaking contributions. 
It is multiplied by an explicit factor of a^, but as seen in the previous term, the power diver- 
gence of the non-oscillatory part of the integral gives rise to an overall scaling of O {g"^)- This 
completes the discussion of the one-loop corrections to the lattice operator for the specific 
displacement vector na used above. It is also straightforward to check the obtained scaling 
of different terms for other choices of the vector na. In general, sub-leading contributions to 
the continuum link are (9(a^), and so by dimensional analysis it has an associated factor of 
momentum squared. On the other hand, it always contains a non-oscillatory term, and as a 
result, the non-continuum contributions and the violations of rotational symmetry scale as 
Oias). 

Given the discussion of the previous paragraphs, we naively conclude that the rotational 
symmetry breaking scales as ~ 0{as) in the continuum limit. It is the one-gluon vertex 
associated with the smeared-operator that is dominating this behavior, with the contribu- 



tions from other diagrams scaling as ~ ag/N for Wilson fermions (eq. (36) and eq. (37)) 
and Us/N"^ from the other loop diagrams compared with ~ from the tree- level match- 

ing. However, this scaling can be further improved by smearing the gauge-field. The 0{as) 
contributions are due to the explicit factor of a? being compensated by a quadratic loop 
divergence, ['K/af', rendering a suppression by only the coupling in the continuum limit, 
analogous to the impact of tadpole diagrams. However, by smearing the gluon field over a 
volume of radius 1/Kg = aNg^ the offending diagrams in fig. [s] scale as 

g 



We have distinguished the smearing radius of the operator, N , from the smearing radius of the gluons, 
Ng, but in principle they could be set equal. 



22 



due to the suppression of the high momentum modes in the gluon propagator. 

The natural question to ask here is what is the scale of the coupling in this process? 
Note that the bare coupling constant of lattice QCD suffers from large renormalization as 
discussed before, so a better-behaved weak coupling expansion of the lattice quantities uses a 
renormalized coupling constant as the expansion parameter. As is suggested by Lepage and 
Mackenzie [501, one first fixes the renormalization scheme by determining the renormalized 
coupling a^g^"' [k*) from a physical quantity such as the heavy quark potential. Then the 
scale of the coupling is set by the typical momentum of the gluon in a given process. In the 
case considered above, the energy scale of the strong coupling constant is dictated by the 
scale of the gluon smearing region as the dominant contribution to the integral comes from 
this region of the integration: k* ~ Tx/{Nga). A better estimate of the scale can be obtained 
by the method explained in Ref. ^50j, but since we are interested in the continuum limit 
where a — ?■ 0, this is already a reliable estimation of the momentum scale of the running 
coupling. 

The analysis in QCD is more complex at one-loop level than in the scalar theory due 
to the presence of the gauge-link required to render the operator gauge-invariant. We have 
found that the contributions from the operator defined in eq. (24) scale in the same way as 
those in the scalar theory, with the violation of rotational symmetry suppressed by factors of 
~ l/N"^, but both tadpole improvement of the extended links and smearing of the gauge-field 
is required. Our analysis of Wilson fermions reveals the contributions to matrix elements 
that violate rotational invariance in the continuum limit at the one-loop level are suppressed 
by factors of ~ cXs/N"^ and ~ ag/N^, and thus for a smearing defined in physical units, 
deviations from rotational invariance scale as (9(a^). Contributions that scale as ~ Us/N 
and are proportional to the Wilson parameter, conserve angular momentum and can be 
absorbed by the operator Z-factor. Most importantly, as in the scalar theory, there are no 
mixings with lower dimension operators that diverge as inverse powers of the lattice spacing. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, a mechanism for the restoration of rotational symmetry in the continuum 
limit of lattice field theories is considered. The essence of this approach is to construct an 
appropriate operator on the cubic lattice which has maximum overlap onto the states with 
definite angular momentum in the continuum. In analogy to the operator smearing proposals 
given in Refs. [23H26] and Refs. [32 | 133 ] , the operator is constructed on multiple lattice sites. 
Using spherical harmonics in the definition of the operator is key to having the leading 
contributions to the classical operator be those with the desired angular momentum. The 
sizes of the contributions are controlled by the scale of the smearing of the operator, with sub- 
leading contributions to both lower and higher dimensional operators that violate rotational 
symmetry being suppressed by - reflective of the pixelation of the operator and fields. 

The A0"^ scalar field theory is shown to preserve this universal scaling of the leading non- 
rotationally invariant contributions at all orders in perturbation theory, compatible with the 
finite size scaling results of A0'^-type theories near their rotational invariant fixed points ^81 
1^ . The same can be shown to be true in gcfr' scalar field theory. 

Gauge invariance somewhat complicates the construction and analysis of analogous op- 
erators in QCD. Although the tree-level lattice operator in QCD exhibits the same scaling 
properties as the scalar operator, extended gauge-links connecting the quark fields generate 
gluonic interactions that contribute to loop diagrams that are power-law divergent. Such 
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10 




FIG. 9. The absolute value of the ratio of the tree-level coefficient, C^q.iq, of a lowest dimension 

(3) 

operator with L = 1 to the tree-level coefficient, Cgg.gQ, of the lowest dimension operator with 
angular momentum, L = 3, resulting from the L = 3 operator in eq. ([T]), as a function of the 
number of included point-shells. 

contributions are either eliminated by tadpole improvement of the extended links, or are 
suppressed by smearing of the gauge-field. We find that it is the physical length-scales and 
continuum renormalization-scale that dictate the size of matrix elements. The leading non- 
continuum corrections from the one-loop diagrams preserve angular momentum, scaling as 
~ agCL for Wilson fermions, and can be absorbed by the operator Z-factor. In contrast, con- 
tributions that violate rotational symmetry are suppressed by a^a^ as a — )■ 0. While we have 
chosen a specific form for the smeared operator, we expect that the results, in particular the 
scaling of the violations to rotational symmetry, are general features of a smeared operator 
with any (smooth) profile. Also, it is worth mentioning that although the calculations pre- 
formed in this work, and the subsequent conclusions, relate operators and matrix elements 
in H(3) to those in 0(3), the methodology and results are expected to hold in relations 
between H(4) and 0(4). Instead of working with operators formed with spherical harmonics 
to recover SO (3) invariance, one would work with operators formed with hyper-spherical 
harmonics to recover 0(4) symmetry. 

We conclude the paper by discussing the practicality of our result for the current LQCD 
calculations as well as its connection to the infra-red (IR) rotational invariance recovery of 
the lattice theories: 

• It is important to understand and to quantify the violation of angular momentum 
conservation in the states and matrix elements calculated using Lattice QCD with the 
lattice spacings currently employed. One interesting result is that by using the tadpole- 
improved operator extended over several lattice sites and built from the smeared gauge 
links, the quantum corrections introduce non-continuum corrections to the tree-level 
results that are suppressed by at least a^, i.e. they do not introduce power-divergent 
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contributions. As an example, suppose that a lattice calculation aims to determine a 
matrix element of an operator with L = 3. Then, as is demonstrated in fig. ([9]), the 
coefficient of the lower dimensional derivative operator with L = 1 is almost 10 times 
larger than the coefficient of the L = 3 derivative operator when the operator is defined 
over one lattice site, = 1. The computational time required to accurately perform 
the subtraction of the L = 1 contribution is significant for a smearing scale of, say, 
A ~ 2 GeV. Fortunately, by halving the lattice spacing and smearing the operator 
over just two point shells (A^ = 2), the contamination from the lower dimensional 
operator is reduced by a factor of ~ 3, requiring a factor of ~ 10 less computational 
resources to accurately perform the subtraction at the same level of precision. Further, 
by smearing the operator over ten point shells, the contamination from the lower 
dimensional operator is reduced to ~ 1% of its value at A^ = 1. Given that the lattice 
spacing associated with A = 2 GeV is a ~ 0.1 fm for A^ = 1, to be able to smear out to 
the N = 2 shell requires a lattice spacing of a ~ 0.05 fm, pushing the limits of current 
lattice generation. To smear out to the A^ = 10 shell would require a lattice spacing 
of a ~ 0.01 fm which is currently impractical. 

• The restoration of rotational invariance as discussed in this paper regards only the 
UV asymptote of the lattice theories: as one reaches a good pixelation of a region of 
space where the lattice operator probes, the identification of eigenstates of the angu- 
lar momentum operator becomes possible. In the other words, the more point-shells 
included in the lattice operator, the larger overlap the operator has onto a definite 
angular momentum state. However, the full recovery of rotational invariance in the 
lattice theories requires the suppression of rotational symmetry breaking contributions 
to the physical quantities not only as a result of short- distance discretization effects, 
but also as a result of boundary effects of the finite cubic lattice in the IR regime of 
the theories. The finite size of the lattice imposes (anti-)periodic boundary conditions 
on the lattice wavefunctions which enforces the lattice momenta to be discretized, 
p = where L is the spatial extent of the lattice and n is a vector of integers. The 
IR rotational invariant theory is achieved as the lattice becomes infinitely large, corre- 
sponding to a large number of point-shells in the momentum space. However, beyond 
this intuitive picture, one needs to examine in a quantitative way how this recovery 
takes place in the large volume limits of the lattice theories in the same way as it was 
discussed for small lattice spacing limit of the theories. One quantitative explanation 
of this IR recovery, has been given recently in Ref. |5T] in the context of the extrac- 
tion of phase shifts in higher partial-waves from the energies of scattering particles 
in a finite volume using Liischers method. The idea is that as one includes higher 
momentum shells, the number of occurrence (multiplicity) of any given irrep of the 
cubic group increases. As a result, for a fixed energy in the large volume limit, linear 
combinations of different states of a given irrep can be formed which can be shown to 
be energy eigenstates; and the energy-shift of each combination due to interactions is 
suppressed in all but one partial- wave in the infinite- volume limit. So, although each 
irrep state has an overlap onto infinitely many angular momentum states, the high 
multiplicity of a given irrep in a large momentum shell generates energy-eigenstates 
which dominantly overlap onto states of definite angular momentum, and the mixing 
with other angular momentum states becomes insignificant in the large volume limit. 
This picture also helps to better understand the mechanism of the UV rotational in- 
variance recovery due to the operator smearing. It is the high multiplicity of the irreps 
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in large (position-space) shells that is responsible for projecting out a definite angular 
momentum eigenstate. These large shells are obtained by reducing the pixelation of 
the lattice by taking a ^ in position-space, or increasing the size of the lattice by 
taking L — > oo in momentum-space - both are required in order to recover rotational 
invariance from calculations performed on a lattice. 
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Appendix A: Operator Basis 

In this appendix, a basis for composite local operators is presented. Any local operator 
that is bilinear in the scalar field with L spatial indices, and that is invariant under cubic 
transformations, can be written as 

= 0Ux) gS^i...,,0(x) , (Ai) 

where Qf^i^,„i^ is a homogeneous function of the operator Vj, and degree d {d > L) is defined 
to be the number of Vs. Their forms are determined by the symmetric traceless tensor of 
rank L that respect cubic symmetry constructed from d Vs. The operators composed of 
fewer than seven derivatives and with no spatial indices are 

(/.t (x)0(x) 
0t (x) V'<f> (x) 
0t(x) (V2)'0(x) 

0t(x) 5]V|0(X) 

i 

(/)t (x) (V2)V(x) 

0t(x)V^ E^J'^W 

j 

4>H^)J2^Ui^) ■ (A2) 
j 

Except for three of these operators which explicitly break the rotational symmetry, they 
transform as L = under rotations. 



0(0) 


(x) 


0(2) 


(x) 




(x) 


q{4,RV) 


(x) 


0(6) 


(x) 


q{6,RV;1) 


(x) 


q{6,RV;2) 


(x) 
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The operators with one spatial index with up to six derivatives are 

0t (x) V.0(x) 
0t (x) V2V,0(x) 
Of (x) = 0t (x) (V2)^v.0(x) 
O)''^^) (x) = 0t (x) J2 V| (x) . 



of 


(x 


of 


(x 


of 


(x 


{5,RV) 


(x 



(A3) 



There is one operator which breaks rotational invariance, and the rest transform as L = 1 
under rotations. 

The operators with two spatial index with up to six derivatives are 



of (x) 



Of (X) 



fx) 



of (x) = 0t (x) (V^) 



2\2 



1 



0(x) 
0(x) 
0(x) 



0(x) 



(A4) 



There is one operator which breaks rotational invariance, and the rest transform as L = 2 

under rotations. 

Operators with three, four and five spatial indices which have L = 3, L = 4 and L = 5 
respectively are listed below. There is no operator which breaks rotational invariance up to 
six derivatives: 



of (x) = (x) 



ViVjVfe - -V^ {SijVk + SjkVi + SkiVj 
o 



0(x) 



of (x) = 0t (x) 



1. 



0(x) , (A5) 



0% (X) 



0% (x) 



it(x) 



ViVjVfeV 

1 



0(x) 



0t (x) V' 



ViVj-VfcVz 



0(x) , (A6) 
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0%m (x) = 0^ (x) [V.V, V,V,V, 



ijklm 

m 

1 2 

+ 7^ (V^) [{5ij5ki + SikSji + Sii5jk) Vm + {SijS 

km ^ik^jm ^im^jk) ^ I 

+ i^ij^ml + ^im^jl + ^il^jm) Vfc + {SimSkl + ^ik^ml + ^il^mk) Vj 
+ {^mjhl + ^mk^jl + (^mi^ifc) Vj]] (x) . (A7) 



Note that as demonstrated in eq. (A2), there can be more than one operator that breaks 
rotational invariance at a given order in derivative expansion. To arrive at a notation that 
is general and useful, one can use the fact that any cubically invariant polynomial of a 
three- vector V, can be expanded in terms of only three cubically invariant structures, 

k k k 

The number of times each structure appears in a derivative operator, as well as the number 
of free indices, uniquely specify the operator. For example, with nine derivatives and one 
spatial index, one can make four independent operators, 

Of^^^ (x) = 0t(x) (V^)V.0(x) 



(^(2,1,0) 



0?'°'^^ fx) 



(x) = 0t(x) {V'Y{y,VI^ V. 



0p'°)(x) = 0t(x)l>JVt| V.0(x) , (A9) 




and generally. 



(x)=(vt(ev')"(e^°) 



o.""'"-" (X) = (v^)- \yyi\ \yyi\ v,4>(^) . (m 



It is then obvious that d = 2m + 4n + 6p + L gives the total number of derivatives in the 
operator, where L is the number of free indices. For n = p = 0, the operator is rotationally 
invariant with angular momentum L. 



Appendix B: Rotational Invariance Violating Coefficients : An Example 

In this appendix, an explicit derivation of a rotational invariance violating coefficient in 
both coordinate-space, and momentum-space formalism, introduced in section |TT| is pre- 
sented. Consider the position space operator ^qq^ (x; a, A^) where superscript indicates that 
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only operators with four derivatives are retained in the expansion of ^oo- The goal is to 
derive the LO correction to the continuum values of coefficients Coo''oo ^^'^ ^ooocT''- 

C (x; a, iV) = (x) [{Na)' cSJoo (V^)' + (Na)' C^lST^ (V^ + Vj + V^)] (x) 

= f ^ E fdyy' f dn, e^^-^-^^t (^) . v)^ (x) Foo (n, 

p " J 

The y integration is 



.(Bl) 



(B2) 



and the coefficients a, /3 and 7 can be determined. It is easy to see that coefficient a makes 
the dominant contribution in the large N limit. Using 



/ [p^) (p.A)^ = f (p^) (p\Af + aJ2 (A^y) 
p p \ i / 



(B3) 



for any rotational invariant function / of the vector p, with 

1 



P 



a 



i^pt - IpI ) , 



(B4) 



one finds that the deviations of Cgo''oo Cqq ocT'* from their continuum values are 



SC, 



(4) 

00,00 



96^/7^ 



Pt^O 



3 COS (27riV|p|) 



SC, 



(4;W) 



00,00 



96^/7^ 



47r2|p|''A^2 
3cos(27riV|p|) 



) (-3pt + IpI^) 

){5pt-\pf) . 



(B5) 



The emergence of rotational invariance violating coefficients from the momentum-space 

construction is somewhat less obvious. From eq. (14) and eq. (15) the operator (k; a, A^) 
can be written as 



ei'J (k; a,N) = <p (k) (-k) {NaT C'Zo Ikp + (NaT C^loo'iK + K + K) 



4^(4;W)/, 4 



p Li,Mi,L2,A^2 

X (LiO; L2O |00) (LiMi; L2M2 |00) Yl,m, (%) >l.m. (^^p) 



X 



dy y^ jLi (a^ |k| y) jL2 {2nN \p\y) 



fc4 



(B6) 



where only the terms of order fc^ are retained from the integral. As such, only Li = 4 with 
y^-i-4 {flp) and i^p), and Li = with Fqo (^p); contribute to the sum. This reduces the 
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relation to 



liJ(k;a,iV) = 6v^ roo(%)>oo(^^p) / dy y 



{aN\k\yy 
120 



■jo(27rAr|p|y) 



+ 9 [(40; 40 |00)' (^fc) {^p) + (40; 40 |00) (44; 4 - 4 |00) ^4 {^t) ^4-4 {ytp 
+ (40; 40 1 00) (4 - 4; 44 100)^-4 (f^^) ^44 iyip)] 



X / 2/ 



2 («A^|k|z/)' 



945 



j4 (27riV|p|y) 



Using the relations 



(ipi'n±4(^^p)) 




16 V 27r 



p p 
and keeping the LO term in 1/A^ from the y integration gives 



fi)(k;a,iV) = 3(aiV|k|)^5^ 

21 A 



cos (27riV|p|) 
47r|p|^A^2 



120 



^00 



+^ - IPl*) 



Finally, we use the relation 

+ + 



3 / 3^ 

16 V vr^^° ^^^^ + 16 V ^ ^^^-^ ^^^^ + ^^^^^ 



^Foo (^^^) + (^^fe) + y ^ (n-4 (^^^) + I44 (fife)) 



(B7) 



(B8) 



(B9) 



(BIO) 



to identify the coefficients ^Cqo^qo ^'^d '^^ooiocT^ from eq. (B9) 

3cos(27riV|p|) 



5C, 



(4) 

00,00 



, /tt .^-^ 



96x/7r 



PT^O 
V p^O 



47r2|p|^iV2 

3cos(27riV|p|) 
47r2|p|^Ar2 



) (-3P^ + IpI^) 
) {^pI - IpH 



(Bll) 



which recovers the position-space results given in eq. (B5). 



Appendix C: Matrix Elements for Non-Zero External Momentum 

The loop calculations presented in the body of this paper have been performed for van- 
ishing external momentum, therefore only the quantum corrections to the L = operator 
have been considered. In this appendix, the generalization to non-zero external momentum 
is presented, where the one-loop correction to the two-point function with an insertion of 
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k+P 



FIG. 10. One-loop contribution to the two-point function with an insertion of the operator in 



the smeared operator is considered in scalar g(j)^ theory, see fig. 10 The loop integral to be 
evaluated is 



LM 



where 



3 d-k 



4vrA^3 



E 



sm 



n a ^ ' 



2 / 



,ik-na 



k"^ + m?] \ {k + P] +m 



k + P 



E- 



sm 



{k^ + Pf,) a 



(CI) 



(C2) 



Note that the operator is smeared over a physical region whose size is small compared to the 
hadronic scale, and as a result the external momenta are small compared to the the scale of 
the operator A = 1/Na. Therefore one may perform a Taylor expansion of the loop integral 
in Pi /A to obtain 



J, 



LM 



3 1 



levr^ A2 



nN 



nN 



dy 3l (qy) 







Ylm (^^c 



-2 



i=l 
3 



m 

\2NJ ' 



xlAN^y sin^ ( ^) + AN' sin^ (^ + 
I \2NJ \2N 2NA 




E 

fc=0 



4iV^ Ell I sin (I) sin (^) + AN' Eli cos (|) sin^ (^) 



4iV2 ^ti sin^ (^) + 4Ar2 sin^ + ^^a) + 



A2 



(C3) 



where q = k/A, ^4 = k^/A, and only the leading term in the Poisson sum is retained. As 
was shown before, the non-zero terms in the Poisson sum are suppressed by at least 1/N' 
compared to the continuum operator insertion in the loop. 

The first term in the above Taylor expansion corresponds to the zero external momentum 
in the loop, therefore at LO, it contributes to the L = operator, and the sub-leading 
rotational invariance breaking operators can be easily shown to be suppressed by 1/N' 



using the procedure described in section |II B[ Note that the loop integrals one needs to deal 
with in g(f)^ are more convergent than comparable integrals in X(j)^ theory, which simplifies 
the discussion of the scaling of the different contributions. 
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The next term in the Taylor expansion of the loop integral can be expanded in large 
since the integral is convergent. The numerator has an expansion of the form 




where the rotational invariance breaking terms are suppressed by at least and the 

leading contribution to the above integral modifies the L = 1 matrix element, while the 
L = term is suppressed by 1/A compared to the L = 1 contribution. The next terms in 
the Taylor expansion give rise to contributions to the L = 2, 3, ... matrix elements at the LO 
in 1/A, while the rotational invariance violating terms remain suppressed by at least 1/A^^ 
compared to the LO contributions. 



Appendix D: Links on the Grid 



In this appendix, the method to evaluate the link at O {g) on a three-dimensional grid is 
outlined through an example, and the result is generalized to other similar cases. The link 
is constructed to be the closest link to the continuum diagonal link in the continuum. 

Suppose that the link lies between points x and x + na on a cubic lattice where: na = 
Oo {Q, 1, 0). Q is an arbitrary integer and oq is a finite number denoting the original lattice 
spacing which is not necessarily small. Then the paths which make minimal area with the 
diagonal path can be formed easily. Among those, the paths which are symmetric under 
reflection around the midpoint of the path are desired since they have somewhat simple 
forms. One such a path in shown in fig. |8^ for Q = 2, where it is straightforward to show 
that: 



U, 



(19) 

(Q,i,o) 



(<?) = igaoe 



iq-na/2 



Ay (g) + (q) 



sin f Qq^ao ^ qyUo 

sin(g^ao/2) 



cos 



+ 



(Dl) 



If the lattice spacing is halved, the closest link to the diagonal path can be obtained by 
adding up two paths each of the form above with an appropriate phase factor and where na 
is replaced by na/2, fig. Isb, 



(Q,i,o) 



in (3^) 



2 sin (^) 



Ay (g) 



sm(g^ao/4) 



(D2) 



This process can be repeated to build extended gauge links on finer grids. For the general 
case, where the original lattice spacing is divided by 2^, it is not hard to show that 



^2^(Q,l,0) 



2^ sin(^) 



Ay (g) 



+ 2A^ (q) cos(^+ 



sin (g,ao/2^+i) 



2A'+2 ' 2-^+1 



. (D3) 
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The continuum limit is obtained by taking K oo, which corresponds to a = 00/2^^ — )■ 0, 
recovering eq. (40). Note that after interchanging the gauge field indices properly, this 
expression is applicable to a class of n vectors with one zero component and ni/nj = Q for 
the ratio of the remaining components. 

The above expression for the gauge link in eq. (D3) can be generalized easily to another 
class of n vectors with one component being equal to Q and the other two components each 
being one. For example for na = oq {Q, 1, 1) one obtains 

sin (^) 



U. 



1,1) 



'^2^' sin(^) 



A, (q) e 



iqyao/2^ 



+ Ay (q) e 



+2A, (q) 



sin (Qg.ao/2^+^) 
sin(g^.ao/2^ 



cos 



Qq X 



+ 



-igzao/2-'' 



2K+2 ' 2^^^+^ 2-^+1 



(D4) 



However, since the vector na is symmetric in its y and z components, the link has to respect 
this symmetry as well. In fact, there exist an equivalent path which arises from the first 
path by interchanging the steps in the y direction and the z direction. Taking an average of 
these two paths gives a link which is symmetric in the y and z components, 

sin {^) 



U. 



(ij?) 

2^(Q,l,l) 



1,1) 



^2^' sin(ii^) 



+ Ay (q) COS 



ao 



+2A, (q) 



sin (Qg.ao/2^+^) 
sin (g,ao/2^+i) 



cos 



Qqxao ^ qytto 



+ 



2K+2 ' 2^+^ 2^+^ 



(D5) 



Taking the K ^ 00 limit of the above link gives rise to the rotational invariant link as well 
as non-continuum corrections which start at O (a^). 

Another class of n vectors are those where two components are equal to Q while the other 
one is equal one. For example for na = ao {Q, Q, 1) the link which is symmetric with respect 
to X and y can be shown to have the form: 



fjiig) 

^2K(QxQxl) 



^,^q.Ax/2!l^i_(^) 



sm 



2(A.(g)cos(||^)+A,(g)cos(||^ 



X cos 



Qqxao Qqyap qzap 
2^+2 ~^ 2-^^+2 2-^+^ 



sin {Qiq^ap + g^ao)/2^+^) 
sin((g^ao + gyao)/2^+i) 

Az (q) 



(D6) 



where U the average of two links which are identical upon interchanging the x and y coor- 
dinate axes. This link recovers the rotational invariant link up to corrections of 0{a^). 

For n vectors with equal components, na = ao{Q,Q,Q), there are six equivalent links 
which are averaged over to obtain 



U. 



(19) 



2K{QxQxQ) 



^ .q.Ax/2!^(^) 



X 



2 Ux{q)oos(^^ 
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\2^+i 2-^+^ 
sin {{Qqxao + Qqytto + Qqzao)/2^'^'^) 
sin {{qxttQ + qyUo + qzao)/2^+'^) 



sin (3^) 

+ Ay (g) cos -—^ + 



+ Az (g) cos 



qxao ^ qyao 



2^ 



2K+1 



COS 



Qqxao _^ Qqyao _^ Qqzao 



2K+2 



2K+2 



2K+2 



(D7) 
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which results in O (a^) corrections to the rotational invariant continuum path. It is the case 
that determining the link for a general extended path is quite involved, but the general trend 
that the deviation from the rotationally invariant continuum path is O(a^) is anticipated. 
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